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Abstract 

Crystal and magnetic structures of a series of novel quantum spin trimer system 
Ca3Cu3_ x Ni x (P04)4 (x=0,l,2) were studied by neutron powder diffraction at the temperatures 
1.5-290 K. The composition with one Ni per trimer (x=l) has a monoclinic structure (space group 
P2 1 /a, no. 14) with the unit cell parameters a = 17.71 A, b = 4.89 A, c = 8.85 A and (3 = 123.84° 
at T=290 K. The (x=2) composition crystallizes in the C2/c space group (no. 15) with the doubled 
unit cell along c-axis. Each trimer is formed by two crystallographic positions: one in the middle 
and the second one at the ends of the trimer. We have found that the middle position is occupied by 
the Cu 2+ , whereas the end positions are equally populated with the Cu 2+ and Ni 2+ for (x=l) while 
in the (x=2) the trimers were found to be of only one type Ni-Cu-Ni. Below T/v = 20 K the (x=2) 
compound shows an antiferromagnetic ordering with propagation vector star {[|, |, 0], [— |, |, 0]}. 
The magnetic structure is very well described with the irreducible representation T2 using both 
arms of the star {k} with the magnetic moments 1.89(1) [is and 0.62(2) \xb for Ni 2+ and Cu 2+ 
ions, respectively. We note that our powder diffraction data cannot be fitted by a model containing 
only one arm of the propagation vector star. The Cu/Ni-spins form both parallel and antiparallel 
configurations in the different trimers, implying substantial effect of the inter-trimer interaction on 
the overall magnetic structure. 

PACS numbers: 75.30.Et, 61.12.Ld, 61.66.-f 
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I. INTRODUCTION 



The low-dimensional magnets have been attracting attention during last years since they 
show new interesting quantum effects and also because they are considered as a model system 
to study very complex phenomena, like high-temperature superconductivity in metal oxides. 
AsCu3(P04)4 (A=Ca, Sr, Pb) is a novel quantum spin trimer system [1] in which the three 
Cu 2+ (S = |) spins are antiferromagnetically coupled giving rise to a doublet ground state, 
as determined by neutron spectroscopy [2]. The trimer clusters form (ID) chains with weak 
but not negligible intertrimer interaction [3, 4] leading to a long-range magnetic ordering 
at T c = 0.91 K, T N = 0.91 K and T N = 1.26 K for A = Ca, Sr and Pb, respectively. 
By substituting a Cu 2+ spin in the trimer by Ni 2+ (S — 1) a singlet ground state could 
be in principle realized offering the observation of the Bose-Einstein condensation (BEC) 
in a quantum spin trimer system similar to the field-induced BEC of the bosonic triplet 
state in the spin dimer system TICUCI3 observed by the inelastic neutron scattering [5]. 
The substitution of Cu by Ni was successfully realized in the Ca 3 Cu3_ x Ni x (P0 4 )4 (x=l,2) 
resulting in the new mixed trimer phase with the structure parameters close to the ones of the 
pristine material (x=0). The magnetic excitations in this series were studied by the inelastic 
neutron scattering [6] making use of the structure data reported in the present paper. The 
observed excitations were associated with transitions between the low-lying electronic states 
of trimers. The nearest-neighbor exchange interactions within the trimers in the (x=l,2) 
compounds were determined to be also antiferromagnetic with J Cu _ Cu = —4.92(6) meV and 
Jcu-m = —0.85(10) meV and an axial single-ion anisotropy parameter Dm = —0.7(1) meV. 
The ground state was found to be doublet, triplet and quintet in the Cu-Cu-Cu, Cu-Cu-Ni, 
and Ni-Cu-Ni trimers, respectively that are the basic constituents of the title compounds. 
The hypothesis of realizing the singlet ground state that motivated the present work was not 
met, but without the detailed structural information the analysis of the magnetic excitations 
could not be performed. In this paper we present the results of the neutron and synchrotron 
x-ray powder diffraction study of the crystal and magnetic structures of Ca3Cu3_ x Ni x (P04)4 



II. SAMPLES. EXPERIMENTAL 



Polycrystalline samples of Ca 3 Cu 3 _ x Ni x (P04) 4 (x=0,l,2) were synthesized by a solid state 
reaction using CuO, NiO, CaC0 3 and NH 4 H 2 P0 4 of a minimum purity of 99.99%. The 
respective amounts of the starting reagents were mixed and heated in alumina crucibles 
very slowly up to 600° C and then annealed at 900° C during at least 100 h, with several 
intermediate grindings. The ac magnetic susceptibility xCO = x'(T) +ix"{T) was measured 
in zero external field with amplitude of the ac field 10 Oe and frequency 1 kHz using Quantum 
Design PPMS station. Neutron powder diffraction experiments were carried out at the 
SINQ spallation source of Paul Scherrer Institute (Switzerland) using the high-resolution 
diffractometer for thermal neutrons HRPT [7] (A = 1.866 A, high intensity mode Ad/d > 
1.8 ■ 10~ 3 ), and the DMC diffractometer [8] situated at a supermirror coated guide for cold 
neutrons at SINQ (A = 4.2 A). All the temperature scans were carried out on heating, x- 
ray synchrotron diffraction measurements at room temperature were done at the Material 
Sciences beam line (MS, SLS/PSI). The refinements of the crystal and magnetic structure 
parameters were carried out with FULLPR0F [9] program, with the use of its internal tables 
for scattering lengths and magnetic form factors. 

III. CRYSTAL STRUCTURE 

Both (x=0) and (x=l) compounds have the same space group P2\/a with the structure 
parameters shown in Table I. The diffraction pattern and the refinement plot for the (x=l) 
sample is shown in Fig. 1. There is a small admixture of two impurity phases: whitlockite 
Ca 19 Cu2(P04)i4 and nickel oxide NiO that are indicated by the additional rows of tics. The 
structure parameters of the whitelockite phase were fixed by the values reported in Ref. [10]. 
The NiO phase has two rows of tics due to additional antiferromagnetic phase. The mass 
fractions of the whitlockite and NiO are 2.2(2)%, 0%, and 3.3(1)%, 0.20(2)% in the (x=0) 
and (x=l) compositions, respectively. An important block of structure is the Cu2-Cul-Cu2 
trimer that is shown in Fig. 2 together with the surrounding P0 4 tetrahedra. The crystal 
structure consists of P0 4 tetrahedra connecting the Cu trimers in chains running along a- 
axis as shown in Fig. 3. Figure 4 shows the view of the trimer connectivity projected roughly 
to the (afe)-plane. 
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One can see that the strongest intertrimer interaction is expected to be along 6-axis giving 
the quasi (ID) trimer chains (...AB...) and (...CD...), where A, B, C and D denote the trimers 
shown in Fig. 4. There are two superexchange paths between the Cul(A) and Cu2(B) ions 
that go through two PCvtetrahedra Cul(A)-0-Pl-0-Cu2(B), and also two similar paths 
between Cu2(A)-Cul(B) ions. In addition, the distance between the Cu2(A)-Cul(B) ions 
is the shortest one (dcu2(A)-Cui(B) — 3.4 A, whereas dcui(A)-Cui(B) — 4.9 A) providing the 
largest dipole interaction. The interaction between the ID chains of trimers along a-axis 
is mediated by the two superexchange paths between Cu2 ions Cu2(A)-0-P2-0-Cu2(C) as 
shown in Fig. 4. The intertimer interaction along c-axis (Fig. 3) is the weakest one since 
the path contains an additional Ca-0 link and the ions are separated by large distance 
dcu2-Cu2 = 6.5 A. 

The Ni-atom in the (x=l) sample occupies the Cu2 positions at the ends of the trimers, 
whereas the middle position Cul is occupied by Cu. The end positions are equally populated 
by both Cu and Ni ions. The occupancy factors can be reliably refined due to significantly 
different coherent scattering lengths of Ni (10.3 fm) and Cu (7.7 fm) nuclei. The fit model 
assumed that we have two Cu atoms and one Ni atom per formula unit allowing them to 
occupy both Cul and Cu2 sites. The refined occupancies are listed in Table I. 

The composition with two Ni atoms (x=2) crystallizes in a different space group C2/c 
with a doubled unit cell along c-axis. The transformation from the P 1 2i/a 1 to the C 1 2/c 1 
structure is given by the matrix A = a, B = b, C = 2c and the origin shift p = a/2. The 
structure solution was done with the FOX program [11] using the synchrotron x-ray diffraction 
pattern collected at the wavelength A = 0.9185 A. Final refinement of the neutron diffraction 
data resulted in the structure parameters listed in Table I. The experimental and the refined 
diffraction patterns are shown in Fig. 5. In spite of the doubled unit cell the density of 
the Bragg peaks is the same as for the (x=0,l) compositions due to C-centered Bravias 
lattice. The mass fractions of the whitlockite and NiO impurities amounted to 6.4(2)% 
and 1.56(3)%, respectively. The crystal structure motif in the (x=2) is very similar to the 
one in the pristine compound and thus all the crystal structure parameters can be directly 
compared (see Table I). The noticeable change in the trimer structure is the decrease in Cul- 
04-Cu2 bond angle. Similar to the (x=l) compound the Ni-atoms predominantly occupy 
the end positions, while the middle position is mainly occupied by Cu. Table II shows the 
average cation-oxygen bond lengths and the distortions of all the polyhedra and the bond 
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valence sum BVS for all cations calculated from the experimental distances using FULLPROF 
[9] suite and the BVS parameters from Ref . [12] . The BVS for the oxygen atoms (not shown 
in the table) are very close to 2. 

IV. MAGNETIC STRUCTURE OF x = 2 COMPOUND 

The magnetic susceptibility data are presented in Fig. 6. The susceptibility of the 
composition with (x=l) that does not exhibit magnetic ordering down to 1.5 K is also shown 
in the plot for comparison. The high-temperature part of the susceptibility (T=75-225 K) 
was fitted to x(T) = C/(T-T C w) + B, where C = N A 2S(S+ l)fi B /3 k is the Curie constant, 
B is a constant background term due to the impurity phases. The fit results are shown in the 
inset of Fig. 6. The paramagnetic spin values per magnetic site Cu 2+ /Ni 2+ calculated from 
the refined values of the Curie constant C amounted to S — 0.92(1) and 0.96(2) that are in 
reasonable agreement with the expected average spin-values per magnetic site 0.67 and 0.83 
for x=l and x=2 compositions, respectively. The broad peak at T N — 20 K is associated 
with a transition to the magnetically ordered state. The high statistics neutron diffraction 
(ND) patterns were collected at the temperatures below (T — 1.5 K) and above (25 K) the 
transition at Tjv using neutron wavelength 1=4.2 A. The low temperature pattern possesses 
many additional magnetic Bragg peaks that prove the presence of long range magnetic 
ordering. The difference pattern ("1.5 K"-"25 K") containing purely magnetic contribution 
is shown in Fig. 7. The temperature scan performed in the temperature range from 1.5 to 
21 K confirms that the magnetic Bragg peaks disappear above 18 K. The Bragg peaks of 
the difference pattern can be excellently indexed in the chemical cell of the x = 2 compound 
with the propagation vector k = [||0], thus proving the antiferromagnetic nature of the 
transition at T^. The powder profile matching refinement (Fig. 7) shows that all the peaks 
are well described in the above model (R wp = 8.9, R exp = 6.7, x 2 — 1-73, XBragg = 1-92). 
The small peculiarities near 29 ~ 33.15°, 81.6°, and 69.4° are due to the slight position 
mismatch of the intense nuclear Bragg peaks for the two temperatures. 
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A. Symmetry analysis 



Using the determined propagation vector we performed the symmetry analysis according 
to Izyumov and Naish [13] to derive possible magnetic configurations for both Ni (8f) and 
Cu (4b) magnetic sites of the space group C12/cl (no. 15). For this purpose we used 
program BASIREP [9] to obtain corresponding basis functions ipj (3m-dimensional vectors) 
in the th unit cell of all atoms of the site (j) with multiplicity m. The magnetic moments 
are obtained by the linear combination of the basis functions: 

s,-o = E c XMl ^ , (1) 

X,k L 

where S^o is a m-dimensional column of spins on the position (j), C\^ L are arbitrary 
mixing coefficients, k L enumerates the arms of the propagation vector star {k}, A runs over 
the basis functions that appear in the decomposition of the magnetic representation. The 
magnetic moments of the atoms displaced by the translation t are obtained by the relation: 

Sj(t) = S j0 exp(27nkt), (2) 

where t is a centering translation (||0) + or chemical cell translations. 

The little group Gk of propagation vector contains two elements {1,1}. The star of 
propagation vector has two arms ki = [||0] and k 2 = [— 1|0] that are related by the 2 y 
symmetry operator. The reciprocal (a*b*)-plane showing the propagation vector star in both 
centered and primitive unit cells is presented in Fig. 8. There exist two one-dimensional real 
irreducible representations T\ and r 2 with the characters (1,1) and (1,-1), respectively. For 
completeness we give the k- vector in the settings used in the Kovalev's book [14]. Kovalev 
uses 5112/6 settings for the space group no. 15 with the transformation matrix: A = a, 
B = c and C = — b, where the capital and lowercase letters are the basis vectors for 
B 1 12/6 and standard C12/cl settings. The relationship between the centered B2/b cell 
and Kovalev's choice of primitive lattice is given by (page 57 in [14]): a = A, b = — (A+C)/2 
and c = (—A + C)/2, where the lowercase letters stand for the primitive cell. Using the 
above matrices we find that the k- vector star in the primitive Kovalev basis reads: {[0, 0, |], 
[0, |,0]}, which corresponds to the star k 4 in the Kovalev's notations. The representation 
of this k 4 -vector group contains two one-dimensional irreps T\ and r 2 (pages 59 and 229 in 
[14]) in accordance with the BASIREP calculations. 
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The decomposition of the axial vector representations for Cu and Ni sites reads: 3r 2 and 
3ti03t 2 , respectively. Only r 2 appears in the decomposition for both magnetic sites and 
hence we conclude that the magnetic ordering goes according to the representation r 2 . In the 
th cell there are two Cu-atoms and four Ni-atoms (centering translation excluded). Both 
Ni and Cu sites are split into two independent orbits as shown in Table III. The atoms on 
the second orbits are obtained from the first ones by applying a rotation 2 y around (0, y,\). 
Table III shows the basis functions of r 2 that will be used below. Actually, in our case 
the basis functions for both orbits and both arms of the star {k} can be chosen to be the 
same as the functions for the orbit 1 and the propagation vector ki. However we have also 
constructed a special case of the basis functions for the orbit 2 shown in Table III for the 
purposes we explain below. 

Inside of each Ni-orbit the spins are antiparallel and have the same magnitude for the 
r 2 representation. The trimers are formed by the two Ni-atoms related by inversion with 
the Cu in the inversion center. Since the Cu atoms are in the positions (0 \ 0) or (0 \ |) 
the inversion about Cu moves the Ni atom out of the th cell to the neighboring cell shifted 
either along x or y-axis (shift along z-axis is not important). This translation reverses the 
Ni spin according to (2) leaving the Ni spins parallel in the trimers. 

B. Magnetic structure determination 

The spin components for both Cu and Ni atoms are unrestricted by symmetry giving in 
general 12 independent parameters: one Ni spin and one Cu spin for each orbit. However, 
we constrain the sizes of the spins of the atoms to be equal in both orbits, because the 
inner-atomic energies generating the atomic spin are much larger than the inter-atomic 
exchange interactions. The atom positions were fixed by the values determined from the 
HRPT diffraction pattern measured at 25 K with A = 1.886 A. We performed a simulated 
annealing minimization [15] of the integrated intensities of the 46 magnetic Bragg peaks 
using FULLPR0F program for this general model (A) using one arm of {k}. Finally the result 
of the simulated annealing search was refined using usual Rietveld refinement of the powder 
diffraction pattern. The best fit (model A) is shown in Fig. 7. The fit quality is not really 
good if we compare the x 2 = 4.5 with the one obtained in the powder matching fit x 2 = 1-7 
as shown in Fig. 7. We have to conclude that there is no satisfactory solution in the model 
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with one arm of the propagation vector star {k}. However, for completeness we present the 
fit results for this model in Table IV. 

We have found a real solution considering both arms of the star {k}, which excellently 
fits to the experimental data (model B in Fig. 7). At first, we considered only orbit 1 with 
the propagation vector ki, i.e. only half of the atoms, and obtained an excellent fit. Then 
we have constricted the basis functions of the vector k 2 from the ones for the vector ki by 
the using the relation [formula (9.15) of Ref. [13]]: 

^(k 2 , f) = exp(-2nik2a p ( y g 2 ,j))Sg 2 R( y h g2 )if)(k 1 J), (3) 

where g 2 is the symmetry element generating the arm k 2 , j and f are the initial atom 
number and the atom number after applying g 2 , a v (g 2 ,j) is a translation returning the 
transformed by g 2 atom j' to the 0^-cell, R{h 92 ) is the rotation matrix of the operator g 2 , 
S 92 is 1/-1 for the proper/improper rotation. The basis functions obtained this way are listed 
in Table III. Using these basis functions and considering only orbit 2 with the same mixing 
coefficients C\^ 2 of formula (1) as for the orbit 1 we get identical Bragg peak intensities. 
Thus, the set of the structure factors is essentially the same for both orbits provided that 
the basis functions are related by the transformation (3) and the propagation vectors by the 
matrix R{h g2 ). We note, that it cannot be a general assertion valid for arbitrary space group 
since a crystallographic site can split into the k-star orbits containing different number of 
atoms. The transformation (3) moves the atoms from orbit 1 to orbit 2 and rotates the Cu 
spin by n around y axis, while for the Ni atoms it reflects the spin about (ac) plane. 

In our case the equivalence of the structure factors for the orbit 1, k x and the orbit 2, k 2 
with the basis functions listed in Table III can be easily seen if we consider the primitive 
unit cell. The primitive basis vectors are related to the ones of the C-centered lattice as: 
a = (A + B)/2, b = (—A + B)/2 and c = C, where the lowercase letters stand for 
the primitive basis. In the Oth primitive unit cell, the spins of Cu and Ni on each orbit 
are parallel and the spins on different orbits are related by a 7r-rotation about B. The 
propagation vectors are k x = [|, 0, 0] and k 2 = [0, |, 0] as shown in Fig. 8. The origin of the 
space group can be shifted by \ along c-axis. Since the structure factor is calculated only for 
the atoms on one orbit we shift orbit 2 by (0, 0, — |), so that the Cu-atom stays in the same 
position (|, |, 0) for both orbits. The Nill and Nil2 are in (x, y, z) and (1 — x, 1 — y, 1 — z); 
the Ni21 and Ni22 are in (1 — y, 1 — x, z) and (y, x, 1 — z). The magnetic structure factor 
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is F(H) oc J2j Mij exp (Hr,), for the scattering vector H = h + k, where h is a reciprocal 
lattice vector of the crystal structure, Mi 3 - = H x [Mj x H]/H 2 and the sum runs over Cu 
and two Ni-atoms at the positions Tj with the magnetic moments Mj. The phase factor for 
Cu is the same for both orbits, the phase factors for Nill and Ni22, and Nil2 and Ni21 will 
be the same if we choose the reflections Hi = (h, k, I) + ki and H2 = (k, h, —I) + k2 for 
the orbit 1 and 2, respectively. The reciprocal vectors Hi and H 2 are related by n- rotation 
about B*-axis (Fig. 8) and hence have the same length and give the Bragg peaks at the same 
2#-position in the powder diffraction pattern. Since the spins on two orbits are also related 
by a 7r-rotation about B-axis, which is collinear to B*, the vectors Mj_j are also related by a 
7r-rotation about B-axis. Hence, the intensity I oc| F | 2 will be the same for the Bragg peaks 
located at Hi and H2 for orbits 1 and 2, respectively, and the powder diffraction patterns 
generated by orbits 1 and 2 will be identical. 

The two orbits do not interfere with each other because of different propagation vectors. 
Hence using all the atoms we naturally get the same fit quality with \[2 smaller mixing 
coefficients. The results of the fit are shown in Fig. 7 (marked as model B). This model 
contains only 6 refinable parameters and gives the same fit quality as the powder matching 
refinement, implying that the fit cannot be better for the given propagation vector k. The 
imperfection of the fits near 26 = 33.2° and 54.7° seen in both the powder match and the 
model B difference curves is apparently due to a "non-ideal" subtraction ("1.5 K-25 K") 
of the large nuclear peaks (002) and (110) at these angular positions. Another explana- 
tion could be the presence of weak ferromagnetism but it is beyond the accuracy of our 
experimental data. 

In this model there is no mixing of the basis functions of ki and k 2 on the same orbit and 
thus it gives constant moment configurations for any direction of the spins. The assumption 
of having the same mixing coefficient for the atoms on the different orbits and belonging to 
the different arms is not dictated by symmetry, because the coefficients C\^ L in formula (1) 
are independent quantities for ki and k 2 . However, in our particular case this assumption 
gives an excellent fit and good spin values as shown in Table IV (Model B). 

The best fit magnetic configuration is shown in Fig. 9. The figure shows 1/4 part of 
the magnetic unit cell. The whole magnetic cell contains 48 magnetic atoms. The mutual 
orientation of the spins in the trimers are different for the trimers on orbits 1 and 2. On 
the 1st orbit the Ni and Cu spins are close to antiparallel configuration, whereas on the 
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2nd orbit they are close to a ferromagnetic coupling. Intuitively one would expect to have 
the same spin orientation in all the trimers, because the intra-trimer interactions should be 
the strongest ones. To make the trimers identical one should constrain the Ni-spin to be 
parallel to c-axis and the Cu-spin to be in the plane perpendicular to the c-axis. For this 
configuration Cu and Ni spins are perpendicular for both trimers. However, this constrained 
model gives a very bad goodness of fit x 2 — 8-67 and has to be rejected. Since the spins in 
the trimers are close to a parallel orientation we also tried to constrain them to be parallel 
(Model C in Table IV). This constrained model gives only slightly worse x 2 than the general 
case. The model contains only 4 adjustable parameters: two spin values and two angles that 
is not bad for describing the intensities of 46 magnetic Bragg peaks. The spins are aligned 
roughly along c-axis with very small canting as shown in Table IV. We note that the sub- 
lattices of the trimers corresponding to the different arms of the star {k} do not interfere 
with each other and the magnetic configuration, in which all the spins in the trimers on the 
same orbit are reversed will give the same Bragg peak intensities. Thus we can have two 
type of domains with the reversed mutual orientation trimer spins on orbit 1 and orbit 2. 

From the above proof of the equivalence of the structure factors for the orbit 1 with ki 
and the orbit 2 with k2 one can see that the model with the atoms from orbit 2 and with 
the propagation vector ki ("orbit2+kl") is not an allowed equivalent solution. We tried 
to find a solution by fitting the data to the "orbit2+kl" model, but the best solution has 
much worse xlragg = 5.1 similar to the Model A. The reason is the difference of the phase 
factors for Ni in the structure factor F(H). For example, the phase factor for Ni-atoms in 
the structure factor for the (|, 0, 0) peak is proportional to sin(7nr) for orbit 1, but to sin(7n/) 
for orbit 2 (we use the primitive cell settings). It is clear that the x and y-coordinates have 
no any symmetry relation with the spin value M±. Thus the intensity of the (1,0,0) peak 
for the orbit 2 and propagation vector k x will depend also on the particular value of y, but 
not only on the M±. 

V. DISCUSSION AND CONCLUSIONS 

The Ni 2+ ion always likes to occupy the end positions of the trimers. The valence of Ni if 
it would occupy the position Cul in the center of the trimer have reasonable value according 
to BVS calculations (Table II), and so in this respect Ni could occupy the middle position. 
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Apparently the total crystal energy is minimized for the Cu being in the middle position. 
It is interesting to note that the distorted square pyramids Cu20 5 and the distorted square 
planes Cul0 4 become much less distorted when going from (x=0) to (x=2) compositions 
(Table II). We note, that the decrease in the distortion of the polyhedron around Cu2-site 
in the x=2 compound where it is occupied solely by Ni-ions is in accordance with the fact 
that Ni 2+ (3c? 8 ) is non-Jahn- Teller active ion. 

The trimers in the (x=l) composition can be of three types Cu-Cu-Cu, Cu-Cu-Ni and 
Ni-Cu-Ni with the statistical populations 25/50/25%, respectively. From the diffraction 
data alone we cannot determine these populations, but from the analysis of the magnetic 
excitations in the trimers [6] the real populations were determined to be 36/28/36%, implying 
that the non-symmetric Cu-Cu-Ni trimer is significantly less populated with respect to the 
"ideal" statistical value. The fact that Ni atom does not occupy the middle position has 
precluded from the realization of the Bose-Einstein condensation in this trimer system, 
however it might be worth trying to make the Ni-substitution in the similar compounds 
A 3 Cu 3 (P0 4 )4 (A=Pb, Sr). 

The antiferromagnetic ordering, which we have found in the (x=2) sample occurs at much 
higher temperature (T/v = 20 K) than in the parent (x=0) compound (TV = 0.9 K). The 
higher T/v might be due to the increase in the dipole interaction strength: Ni-ion has two 
times larger spin value and the intertrimer distances between the Cu2-sites along a-axis and 
^-axis are decreased from 4.8 A to 4.5 A, and from 6.5 A to 6.4 A, respectively. In addition, 
the superexchange (SE) coupling is expected to be larger in (x=2) compound. The SE 
interaction between the (ID) trimer chains along 6-axis (as explained in Sec III) is mediated 
by the Cu2-0-P2-0-Cu2 path. The completely closed P 5+ 2p-shells provide SE path that 
can be both anti- and ferromagnetic. The average Cu-O-P bond angle is increased from 
126.3° to 128.5°, implying that the antiferromagnetic SE via Ni(3d)-0(2p) orbitals can be 
larger in the (x=2) compound. 

The antiferromagnetic structure can be well described only using both arms of the prop- 
agation vector star. It is quite unusual that the two-k case can be revealed from the unpo- 
larized powder diffraction data analysis. However, the two-k solution excellently describes 
the data with minimal number of the refined parameters. For further verification of our 
magnetic structure model the single crystal diffraction experiments might be useful. Ac- 
cording to our model the spin orientation in the Ni-Cu-Ni trimers can be both anti- and 
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ferromagnetic. The chains of the trimers running along the a-axis are of two types, one 
consisting of the antiferromagnetic (AFM) trimers and another one with the ferromagnetic 
(FM) trimers. The Ni 2+ ordered magnetic moment 1.9 /j,b is close to the saturated value, 
whereas the Cu 2+ moment 0.6 fis is substantially smaller than the spin-only value. This 
might be due to frustration of the Cu-moment, i.e. some trimers in the AFM trimer chain 
have Cu-spins aligned ferromagnetically with Ni ones and vice versa for the FM-chains. 

In conclusion, we have successfully synthesized and studied the crystal and magnetic 
structures of the novel mixed spin trimers Ca 3 Cu3_ x Ni x (P0 4 )4 (x=0,l,2) by means of neu- 
tron diffraction in the temperature range 1.5-290 K. Our work forms an important ground 
for the inelastic neutron scattering study of the dynamic magnetic properties of this system. 
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TABLE I: The structure parameters and Cu(Ni)-0 interatomic distances in CasCu3(P04)4 (x=0), 
and Ca 3 Cu 2 Ni(P04)4 (x=l) [sp.gr. Pl2 1 /al (no. 14)] and and Ca 3 CuNi2(P0 4 )4 (x=2) [sp.gr. 
C12/cl (no. 15)]. The WyckofT positions are (2a) for Cul, (2b) for Cal and (4e) for other 
atoms for the compounds with (x=0,l) and (4b) for Cul, (4e) for Cal and (8f) for other atoms 
for the (x=2) one. The data are refined from the powder neutron diffraction patterns measured at 
HRPT/SINQ with wavelength A = 1.886 A. Bragg reliability factor -RBragg for the main phase and 
conventional reliability factors for the whole pattern R wp , R exp ,x 2 are also given. The notation of 
the oxygen atoms around the Cu2-Cul-Cu2 trimer are given in Fig. 2. The bond lengths are given 
in A, the angles in degrees and the isotropic thermal displacement parameter B in A 2 . For the 
Cul and Cu2 positions the refined occupancies Cu/Ni are given. 





x=0 


x=l 


x=2 


a, A 


17.62154(8) 


17.71388(9) 


17.7174(1) 


b, A 


4.90205(2) 


4.88512(2) 


4.82109(4) 


c, A 


8.92224(5) 


8.84635(5) 


17.8475(1) 


7, deg 


124.0744(3) 


123.8436(3) 


123.6373(5) 


V, A 3 


638.39 


635.81 


1269.22 


x,y,z 


0,0,0 


0,0,0 


|,o,o 


B (Cul) 


0.64(4) 


0.74(6) 


0.81(8) 


Cu/Ni 


1/0 


0.980/0.020(15) 


0.84/0.16(2) 


x,y,z 


0.1198(1),0.4771(3),0.9430(2) 


0.1213(1),0.4717(3),0.9461(2) 


0.6199(1),0.4638(4),0.4671(1) 


B (Cu2) 


0.39(4) 


0.70(3) 


0.67(4) 


Cu/Ni 


1/0 


0.510/0.490(8) 


0.08/0.92(1) 


x,y,z 


2 ' u ' 2 


2 ' u ' 2 


0,0.951(l),i 


B (Cal) 


0.44(6) 


0.71(8) 


0.42(9) 


x,y,z 


0.2659(2),0.4630(6),0.7265(3) 


0.2648(2),0.4648(6),0.7267(4) 


0.7667(2),0.4441(8),0.3614(2) 


B (Ca2) 


0.80(5) 


0.63(6) 


0.69(7) 


x,y,z 


0.5934(1),0.9897(5),0.2487(3) 


0.5930(2),0.9857(7),0.2480(3) 


0.0933(2),0.9915(7),0.1257(2) 


B (PI) 


0.39(4) 


0.67(4) 


0.22(5) 


x,y,z 


0.8406(2),0.0161(5),0.2164(3) 


0.8417(2),0.0175(5),0.2194(3) 


0.3406(2),0.9498(7),0.1111(2) 
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B (P2) 


0.40(4) 


0.36(5) 


0.35(6) 


x,y,z 
B(01) 


0.6794(2),0.9046(4),0.4240(3) 
0.68(4) 


0.6798(2),0.9023(5),0.4267(3) 
0.68(4) 


0.1793(2),0.0811(6),0.2121(2) 
0.35(5) 


x,y,z 
B (02) 


0.0087(1),0.5519(5),0.2534(3) 
0.68(4) 


0.0098(1),0.5524(5),0.2536(3) 
0.61(6) 


0.5111(2),0.5214(7),0.1301(2) 
0.99(6) 


x,y,z 
B (03) 


0.6006(1),0.2866(4),0.2004(3) 
0.73(4) 


0.5997(2),0.2810(5),0.1927(4) 
0.85(4) 


0.1017(2),0.6968(6),0.0995(2) 
0.68(6) 


x,y,z 
B (04) 


0.5812(1),0.8074(4),0.0918(3) 
0.44(4) 


0.5846(2),0.7987(5),0.0962(4) 
0.76(5) 


0.0801(2),0.1812(6),0.0477(2) 
0.45(6) 


5 (05) 


0.8978(2),0.8549(4),0.3861(3) 
1.09(4) 


0.9014(2),0.8629(5),0.3913(4) 
1.62(6) 


0.3980(2),0.1037(5),0.1972(2) 
2.04(8) 


x,y, z 
B (06) 


0.8542(1),0.9282(4),0.0661(3) 
0.66(5) 


0.8518(2),0.9277(5),0.0641(3) 
1.07(5) 


0.3563(2),0.0482(5),0.0381(2) 
1.11(6) 


x,y, z 
B (07) 


0.3574(1),0.1752(4),0.2467(3) 
0.82(4) 


0.3586(2),0.1658(5),0.2433(3) 
1.20(6) 


0.8542(2),0.1343(6),0.1220(2) 
1.01(6) 


B (08) 


0.7396(1),0.9661(5),0.1521(3) 
0.83(4) 


0.7429(2),0.9709(6),0.1622(3) 
0.95(5) 


0.2412(2),0.9966(6),0.0806(2) 
1.43(6) 


-^Bragg, % 
.R w p , Rexp t X 


1.84 

2.48, 0.92, 7.23 


1.68 

2.68, 1.11, 5.82 


1.69 

3.0 , 1.57, 3.68 


Cul-04 
Cu2-04 
Cul-Cu2 
Cul-04-Cu2 


1.917(2) 
2.086(4) 
3.534(2) 
123.9(2) 


1.918(2) 
2.099(4) 
3.556(2) 
124.5(2) 


1.939(3) 
2.045(4) 
3.352(2) 
114.6(2) 


Cul-03 
Cu2-02 
Cu2-07 
Cu2-08 
Cu2-06 


1.970(2) 
1.942(2) 
2.175(4) 
1.898(2) 
2.051(3) 


1.953(2) 
1.989(2) 
2.023(3) 
2.121(4) 
1.946(2) 


1.938(2) 
2.003(3) 
2.055(5) 
1.984(3) 
2.060(4) 
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TABLE II: The average cation-oxygen bond lengths d, the rms distortion of the polyhedra Sd/d 
in units 10~ 4 and the bond valence sum BVS calculated from the experimental distances using 
FULLPRDF [9] suite and the BVS parameters from [12]. C and V are the coordination of the 
polyhedra and the nominal cation valence, respectively. 





x = 


x = 1 


x = 2 


At. 


C 


V 


d, A 


Sd/d 


BVS 


d, A 


Sd/d 


BVS 


d, A 


Sd/d 


BVS 


Cul 


4 


2 


1.9450(9) 


1.5 


1.953(5) 


1.9358(10) 


0.5 


2.000(5) 


1.9387(13) 


0.002 


1.983(8) 


Nil 


4 


2 












1.869(5) 






1.853(41) 


Cu2 


5 


2 


2.0290(13) 


24.3 


2.013(6) 


2.0336(13) 


10.9 


1.950(12) 






1.947(41) 


Ni2 


5 


2 












1.822(12) 


2.0293(16) 


2.2 


1.819(8) 


Cal 


6 


2 


2.3331(9) 


5.1 


2.253(6) 


2.3273(8) 


4.4 


2.286(5) 


2.3361(16) 


18.4 


2.290(11) 


Ca2 


3 


2 


2.3714(21) 


2.0 


1.010(6) 


2.3707(22) 


6.0 


1.020(6) 


2.4017(25) 


7.4 


1.255(8) 


PI 


4 


5 


1.5384(17) 


2.7 


4.958(22) 


1.5372(18) 


1.3 


4.969(23) 


1.5298(23) 


2.7 


5.074(31) 


P2 


4 


5 


1.5343(17) 


2.5 


5.014(22) 


1.5357(18) 


4.9 


5.005(24) 


1.5297(23) 


2.6 


5.077(32) 
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TABLE III: Positions of the magnetic atoms in the unit cell, symmetry operators of Gk and 
G (space group G = C12/cl, propagation vector k = [|?|0]) and basis functions for irreducible 
representation r 2 . The Ni-atom is in general position (8f) with coordinates (0.62065, 0.53530, 
0.96795), the Cu-atom is in(4b) position (0, 5, 0), The atoms of the orbit 2 are generated from 
the respective atoms of orbit 1 by the symmetry element (2 0, y, \). The basis functions for orbit 
2 are constructed using formula (3). 



Atom 


eqv. pos. 


sym. op. G;G k 


T2 




T " 

T 2 


orbit 1, ki = [||0] 










Mil 


x,y,z 


i;i 


100 


010 


001 


Nil2 


-x+l,-y+l,-z+l 


i;i 


-100 


0-10 


00-1 


Cul 


0,§,0 


i;i 


100 


010 


001 


orbit 2 (2 0,y,|), k 2 = 


[-Ho] 








Ni21 


-x+l,y,-z+| 


2 0,y,|; 1 


100 


0-10 


001 


Ni22 


x-i,-y+i,z-| 


c x,0,z; -1 


-100 


010 


00-1 


Cu2 


u >2 > 2 


2 0,y,|; 1 


-100 


010 


00-1 
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TABLE IV: Magnetic model parameters for Ca3Cu3_ x Ni x (P04)4(x=2) refined from the diffraction 
data shown in Fig. 7. The numeration of the atoms is the same as in Table III. M is the size 
of the magnetic moment, (j) and 9 are spherical angles with c (azimuth) and b (zenith) axes, 
respectively. The graphical illustration of the spherical angles is given in Fig. 9. The errorbars are 
given only for the independently refined parameters. The model A does not fit the data, but given 
for completeness. See the text for details. In the models B and C (f>Ni2i = <Azvui) 6ni21 = ^ — Onui, 
4>Cu2 = 4>Cu\ + 9cu2 = Qcui- In the model C the spins of Ni and Cu are constrained to be 
(anti) parallel in the trimers. 





model A 


model B 


model C 




M, hb 




9 


M, flB 




9 


M, fi B 


4> 


9 


Mil 


1.760(2) 


104.3(1.2) 


78.4(1.5) 


1.892(9) 


176.1(8) 


83.98(45) 


1.883(8) 


175.9(9) 


83.4(4) 


Cul 


1.196(4) 


21.2(4.4) 


47.5(3.1) 


0.62(2) 


153.5(3.7) 


103.3(2.9) 


0.57(1) 


175.9 


96.6 


Ni21 


1.760 


12.4(1.7) 


69.9(1.3) 


1.892 


176.1 


96.02 


1.883 


175.9 


96.6 


Cu2 


1.196 


26.9(4.3) 


127.7(4.1) 


0.62 


333.5 


103.3 


0.565 


355.9 


96.6 






14.1, 6.7, 4.4, 5.1 


8.9, 6.7, 


1.76, 1.93 




9.2, 6.7, 


1.86, 2.05 
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FIG. 1: The Rietveld refinement pattern and difference plot of the neutron diffraction data for 
the sample Ca3Cu3_ x Ni x (P04)4 (x=l) at T=290K measured at HRPT with the wavelength A = 
1.886 A. The rows of tics show the Bragg peak positions for the main phase, and two impurity 
phases: whitelockite, NiO nuclear and NiO magnetic peaks (from top to bottom). 
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FIG. 2: (Color online) Fragment of the crystal structure showing Cu2-04-Cul-04-Cu2 trimer 
and the surrounding PO4 tetrahedra. The positions of the atoms and some selected interatomic 
distances are listed in Table I. 
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FIG. 3: (Color online) Projection of the crystal structure of Ca3Cu3„ x Ni x (PO"4)4 on the (ac) plane 
showing connectivity of the trimers. The Cul and Cu2 positions are shown by blue and white 
balls, the tetrahedra are PO4, the red sticks indicate Cu-0 bonds, the big green balls are the Ca 
atoms. 
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FIG. 4: (Color online) View of four trimers (A, B, C and D) projected approximately to the (ab) 
plane in Ca3Cu3_ x Ni x (P04)4 showing possible superexchange paths between the trimers in [ab) 
plane. The Cul and Cu2 positions are shown by blue and white balls, the tetrahedra with yellow 
balls in the center are PO4, the red sticks indicate Cu-0 bonds, the Ca atoms are not shown. Each 
trimer is formed by the central Cul atom and two Cu2 atoms related by inversion with respect to 
the Cul position. 
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FIG. 5: The Rietveld refinement pattern and difference plot of the neutron diffraction data for 
the sample Ca3Cu3_ x Ni x (P04)4 (x=2) at T=290K measured at HRPT with the wavelength A = 
1.886 A. The rows of tics show the Bragg peak positions for the main phase and two impurity 
phases. 
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Ca3Cu,Ni x (P0 4 ) 4 




50 100 150 200 250 300 

Temperature, K 

FIG. 6: (Color online) Real part of the ac magnetic susceptibility is shown as a function of tem- 
perature for the samples Ca3Cu3_ x Ni x (P04)4 with x=l and 2. The inset shows the inverse sus- 
ceptibility with the least-square fit to the Curie- Weiss law. The susceptibility is given per mole 
of the magnetic ions (Ni or Cu). The refined Curie- Weiss transition temperatures Tew an d the 
Curie constants C are show in inset. 
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FIG. 7: The Rietveld refinement pattern and difference plot of the difference magnetic neutron 
diffraction pattern between 1.5K and 25K for the sample with Ca3Cu3_ x Ni x (P04)4 (x=0.2) at 
T=290K at DMC with the wavelength A = 4.2 A. The difference curves are shown for 3 differ- 
ent refinements: profile matching mode (the model contains only unit cell parameters and the 
propagation vector) and two different models. See the text for details. 
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FIG. 8: The reciprocal (a*b*)-plane showing both conventional (direct space centered C12/cl 
setting) and primitive reciprocal cells by the dotted and solid lines, respectively. The primitive 
basis vectors are related to the conventional ones as a* = A* + B*, b* = —A* + B*. The 
propagation vector star is shown by {ki, k^}. 
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FIG. 9: The Oth unit cell of Ca3Cu3- x Ni x (P04)4 (x=2) showing the configuration of the Ni and Cu 
spins. Some of the Ni-spins from the neighboring cells are shown for better visibility of the trimers. 
The unit cell constants at T=25 K are a = 17.724 A, b = 4.815 A, c = 17.836 A, /3 = 123.756° 
(C2/c space group). The spins in the middle of the cell along c direction between the dotted 
lines belong to the orbit 2 (Cu2, Ni21 and Ni22 spins) and have propagation vector k2 = [— ^0], 
the other spins belong to the orbit 1 (Cul, Nill and Nil2 spins) and have propagation vector 
ki = [|§0]. The structure corresponds to the model B shown in Table IV. The crystal axes and 
the spherical angles used in Table IV are also shown. 
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